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Introduction 

The  long-pulse  induction  linac  'aTNRL  [l-sfis  considered  for  the  first  stage  of  a 
two-stage  free  electron  laser  (FEL)  oscillator.  This  research  effort  studies  several 
effects  that  can  significantly  influence  the  performance  of  such  an  FEL.  The  NRL  FEL  is 
parameterized  in  a  set  of  dimensionless  variables  that  can  summarize  several  physical 
effects  without  the  use  of  detailed  calculation  as  well  as  relating  the  physical  to  other 
FEL  designs.  A  waveguide  analysis  shows  the  primary  modifications  on  the  FEL 
interaction,  electron  beam  distribution  functions  representing  energy  spread  and 

emittance  are  evaluated  in  the  high  gain  regime,  and  a  multimode  analysis  of  the 

c 

trapped-particle  instability  is  performed  for  parameters  describing  the  NflL  FEL.  The 
research  in  intended  to  extend  the  simulation  theory  of  high-gain  FEL  oscillators. 

In  the  two-stage  FEL,  the  first-stage  uses  a  normal  FEL  interaction  with  the  usual 
static  undulator  to  produce  an  intense  electromagnetic  wave  in  a  high-power  resonator. 
The  wavelength  in  the  first  stage  is  around  1cm,  and  provides  the  periodic  undulator  in 
the  second  stage  of  the  interaction.  The  electromagnetic  undulator  has  a  shorter 
wavelength  with  a  more  intense  field  than  can  be  attained  from  a  static  magnetic  field, 
and  allows  the  second  stage  to  reach  short  optical  wavelengths  with  a  low  energy 
electron  beam.  But,  the  performance  of  the  two-stage  FEL  is  completely  determined  by 
the  all  important  first-stage.  For  this  reason,  the  emphasis  of  this  research  is 
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concentrated  on  understanding  the  high  power  operation  of  the  high-gain  FEL  oscillator. 

Data  and  design  parameters  from  the  first  six  references  below  provide  a  physical 
description  of  the  NRL  FEL.  Many  physical  variables  like  current  and  beam  energy 
changed  from  experiment  to  experiment  but  remained  in  the  general  range  described 
here.  The  NRL  FEL  uses  a  ymc2  =  7 keV  electron  beam  so  that  y=2.4.  The  beam 
current  /  is  about  200A  in  a  beam  radius  of  re  =  0.5cm.  Reference  [4]  gives  the  angular 
spread  of  the  electron  beam  as  A8  =  0.014,  and  there  is  a  2%  energy  spread.  The 
undulator  has  a  X0  =  4cm  period  over  N  =22  periods,  and  the  peak  field  strength  is 
B  =  625Gauss  with  helical  polarization.  The  "undulator  constant"  for  such  an  undulator 
is  K  =  eB\0/2nmc2  =  0.23.  An  a  =3cm  diameter  waveguide  confines  the  radiation  to  a 
small  volume  around  the  co-propagating  electron  beam.  The  lowest  order  mode  has  a 
beam  radius  of  w0  =  0.8cm.  The  combined  waveguide,  resonator,  and  output  coupling 
loss  each  pass  n  is  about  30%,  so  that  in  the  absence  of  gain  the  optical  field  would 
decay  as  ~  e~nl where  Q  =3-»4.  The  long-pulse  NRL  accelerator  provides  current 
for  about  2ps,  or  about  n  =  200  passes. 

The  filling  factor  F  =nr2/KWQ  =0.4  describes  the  reduced  coupling  due  to  the 
mismatch  in  the  electron  beam  and  radiation  beam  sizes.  The  Rayleigh  length  for  the 
radiation  that  would  be  emitted  from  the  electron  beam  with  no  waveguide  is  much 
shorter  than  the  interaction  length  L  =  ArAt)  =  88cm.  It  is  shown  below  that  this  could 
significantly  reduce  coupling  in  the  NRL  FEL.  The  general  criteria  for  the  use  of  a 
waveguide  in  an  FEL  is  derived. 

The  electron  beam  density  p  =  5x1010cm-3  and  the  filling  factor  determine  the 
dimensionless  beam  current  j  =  8N(enKL)2pF/y3mc2  =  320.  When  j  »  1,  the  FEL  has 
the  potential  for  high  exponential  gain  if  there  is  sufficient  beam  quality.  The  weak  field 
power  gain  is  given  by  G  =e^l2')  ^/9  and  is  G  =  1 03,  or  G  =30dB.  These  levels  of 
amplification  are  in  rough  agreement  with  the  experiment,  but  there  is  some  degradation 
due  to  beam  quality.  The  gain  spectrum  bandwidth  is  given  by  A0  =4 j  =  2k  for  this 
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FEL. 

For  a  low-energy,  high-current  electron  beam,  interparticle  Coulomb  forces  can 
cause  plasma  oscillations  to  modify  the  simple  FEL  interaction  mechanism.  The 
relativistic  plasma  frequency  is  0.p  =  (Ane2plyzm)v2  =  3x1 09s-1  .  The  number  of  plasma 
oscillations  during  a  passage  of  the  beam  through  the  undulator  in  the  interaction  time 
L/c  is  A/plasma  =CIpL/2kc  =1.6.  This  is  larger  than  for  most  FEL  experiments,  but  still 
small  enough  that  the  main  interaction  would  be  determined  by  the  exponential  growth 
described  above. 

Another  unique  feature  of  the  NRL  experiment  is  that  it  did  not  use  a  uniform  guide 
field  to  contain  the  low-energy  electron  beam  along  the  undulator.  The  helical  undulator 
field  focuses  in  both  transverse  dimensions  and  causes  N^  =  NK /y  =  2.2  betatron 
oscillations  along  the  interaction  length.  Both  the  radius  of  the  electron  beam  and  the 
angular  spread  contribute  to  the  spread  in  electron  phase  velocities  during  the 
interaction.  Energy  spread  also  deters  bunching.  To  compare  these  effects  on  the 
same  footing,  they  are  all  related  to  the  gain  spectrum  bandwidth  M2N .  For  an  angular 
spread  of  A0  in  a  beam  of  radius  re ,  the  spread  in  electron  phase  velocities  is  given  by 

Op  =  4rc/V  (K2k§r2  +  y2Q2  )/('\+K2)  . 

When  Op  =  7c,  the  spread  in  phase  velocities  is  comparable  to  the  gain  spectrum 
bandwidth,  and  gain  degradation  becomes  important.  Both  the  angular  spread  and 
beam  size  contributions  should  be  matched  consistent  with  the  beam  emittance  in  r.per 
to  optimize  performance.  In  the  NRL  FEL,  cp  =  8  indicating  there  will  be  some  loss  of 
gain  due  to  beam  quality.  More  importantly,  it  is  reported  in  several  of  the  references 
that  the  position  of  the  electron  centered  can  be  off  of  the  undulator  axK  by  as  much  as 
the  beam  radius  0.5cm.  This  causes  a  shift  and  a  spread  in  the  phase  velocity 
distribution  of  the  .beam.  The  shift  alone  is  Aop  =  35  and  should  cause  dramatic 
changes  in  the  interaction  strength.  In  the  high-gain  regime,  the  shape  of  the  electron 


-4- 


distribution  function  becomes  extremely  important  in  determining  the  interaction 
characteristics  along  the  undulator.  A  much  more  detailed  analysis  of  these  features  is 
covered  in  an  appendix. 

FELs  using  Waveguides 

When  an  FEL  uses  a  waveguide  to  confine  the  radiation  along  the  undulator,  there 
are  changes  to  the  simple  FEL  interaction.  Usually  the  goal  of  the  design  is  to  make  the 
changes  as  small  as  possible  while  confining  the  wave  around  the  co-propagating 
electron  beam. 

General  Criteria  for  the  FEL  Waveguide:  When  the  FEL  optical  wavelength  is  long,  free 
space  diffraction  can  spread  the  optical  wavefront  away  from  the  co-propagating 
electron  beam  and  reduce  coupling.  The  natural  distance  for  a  light  beam  of  wavelength 
X  to  double  its  initial  beam  area,  twq  ,  is  the  Rayleigh  length  Zq  =  kwq/X.  Comparing 
this  length  to  the  undulator  length  L  defines  the  dimensionless  Rayleigh  length 
z0  =  kwq  /LX.  Without  a  waveguide,  the  optical  wavelength  is  determined  by  the 
resonance  condition,  X  =  At)(1  +K2)/2y2  where  K=eBX0/2twic2  is  the  undulator 
parameter,  X0  =  kcj2n  is  the  undulator  wavelength,  and  ymc 2  is  the  electron  beam 
energy.  The  Rayleigh  length  becomes  z0  =  2rov02Y2/(1-bK'2)/VX^,  and  expresses  how  an 
FEL  with  a  low  energy  electron  beam  gives  a  short  Rayleigh  length. 

A  design  attempting  to  compensate  for  low  y  by  expanding  the  optical  beam  size 
w 0  would  suffer  reduced  coupling  because  of  a  small  filling  factor  F0  =  nr2iwQ  ;  re  is  the 
electron  beam  radius.  When  natural  diffraction  spreads  the  optical  wave  away  the 
electron  beam  the  filling  factor  is  reduced  further;  assuming  the  light  remains  in  the 
lowest  order  Gaussian  mode  with  its  waist  at  the  center  of  the  undulator,  the  average 
filling  factor  over  the  undulator  length  is  F  =  /r0/(1+1/12zo  ),  or  F  =  o2/(z0+1/12z0)  where 
°e  =re(nJLX)v2  is  the  dimensionless  electron  beam  radius.  Typically,  ae  =  1  in  an  FEL, 
but  the  following  arguments  are  independent  of  that  value  and  only  depend  on  the 
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Rayleigh  length  z0.  For  z0->0,  F  -»0  and  coupling  -»0,  because  the  optical  wave 
spreads  away  from  the  electron  beam  at  the  two  ends  of  the  undulator;  for  z0-m> o, 
F  ->  0  and  coupling  -»0  again,  because  the  wavefront  is  too  large  compared  to  electron 
beam.  The  maximum  filling  factor,  Fmax=F0l 2,  or  Fmax=^3a2,  occurs  at 
z0max  =  (12)_1/2.  The  maximum  is  relatively  broad  in  z0  and  F  max  drops  to  half  of  its 
peak  value,  ^3o2/2,  at  z0  =  0.1 1  and  0.75. 

The  filling  factor  could  be  increased  for  an  optical  wavefront  of  any  size,  if  the 
electron  beam  size  is  increased.  But,  as  the  electrons  move  off  of  the  undulator  axis, 
approaching  the  undulator  magnets,  the  field  strength  B  increases  and  deflects 
electrons  back  towards  the  axis  causing  betatron  oscillations.  The  extra  transverse 
motion  decreases  z  velocity  of  the  relativistic  electrons  and  changes  the  resonance 
condition.  The  change  in  the  electron  phase  velocity  in  a  matched  beam  with  radius  rt 
is  Avp  =  4nN (Kk0ref/(\+Kz).  Gain  degradation  begins  when  the  beam  radius  is  large 
enough  to  cause  Avp  =  7t.  The  limit  on  the  beam  radius  for  good  coupling  is  then 
re  <  ( 1 +K 2) 1  ^Xq/AkKN  1/2.  Generally,  a  typical  FEL  has  N  =  102  so  that  the  electron 
radius  is  restricted  to  a  much  smaller  value  than  the  undulator  wavelength.  For  K  =  1 , 
we  have  re  <  0.01  Xq,  or  for  K  «  1 ,  we  have  re  <  0.01  tyK . 

The  combination  of  (1)  the  restriction  on  the  spread  of  electron  phase  velocities,  or 
beam  radius,  and  (2)  the  requirement  of  good  coupling,  or  large  filling  factor  F  with  a 
given  optical  beam  size  w0,  gives  a  restriction  on  the  Rayleigh  length  in  terms  of  y  and 
N.  To  relate  re  and  h>0,  assume  that  the  filling  factor  is  not  too  small,  say  F  >  F0  >  0.1 , 
so  that  w0  <  3 re .  Then,  the  dimensionless  Rayleigh  length  is  limited  by  z0  <  y2/nK2N2. 
Either  a  low  energy  beam  or  a  long  undulator  can  limit  z0  to  a  small  value  and  decrease 
the  filling  factor.  The  limit  is  relaxed  when  K  is  small,  because  the  electron  beam  can  be 
expanded  to  support  a  wide  optical  wavefront  without  much  diffraction.  Gain 
degradation  begins  to  occur  when  the  coupling  is  reduced  by  small  values  of  F .  The 
peak  value  of  F ,  ^3o2,  decreases  to  less  than  half  its  value  when  the  limit  above 
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restricts  z0  <  0.1 .  Therefore,  an  FEL  design  requires 

Y<  KN/2 

so  that  natural  diffraction  does  not  significantly  decrease  the  interaction  strength.  An 
FEL  with  low  electron  beam  energy  and  a  long,  strong  undulator,  requires  a  waveguide 
for  good  coupling  strength. 

FEL  Resonance  in  a  Waveguide:  When  a  waveguide  is  used  to  confine  the  radiation  near 
the  co-propagating  electron  beam  there  is  a  change  in  the  FEL  resonance  condition.  In 
a  single  waveguide  mode,  the  cross-section  of  the  mode  can  be  used  in  the  filling  factor 
to  estimate  coupling  to  the  mode.  The  electron  beam  size  should  be  close  to,  but 
smaller  than,  the  radiation  cross-section  for  best  coupling.  Assuming  the  electron  beam 
is  on-axis,  the  mode  should  have  a  transverse  electric  field  on-axis  as  well.  Higher  order 
modes  will  average  to  smaller  coupling  if  the  electron  beam  size  is  not  much  smaller 
than  the  mode.  At  NRL,  the  waveguide  is  cylindrical  in  cross-section  with  a  a  =  1.5cm 
radius.  Since  the  electron  beam  size  is  about  re  =  0.5cm,  and  is  comparable  to  the 
fundamental  mode  size,  no  more  than  2  to  3  nodes  should  not  couple  well. 

The  time  dependence  of  the  waveguide  fields  is  taken  to  be  with 

-f  ilry 

longitudinal  dependence  «  e~  .  Waveguide  modes  are  separated  into  two  classes:  TE 
-  transverse  electric  modes  where  the  longitudinal  component  of  the  electric  field  El-  0 
everywhere,  and  TM  -  transverse  magnetic  modes  where  the  longitudinal  component  of 
the  magnetic  field  B2  =  0  everywhere.  At  the  waveguide  wall,  the  TE  mode  boundary 
condition  is  Bz '  =  0  while  the  TM  mode  boundary  condition  is  Ez  =  0.  The  waveguide 
cross-section  and  boundary  conditions  specify  and  eigenvalue  problem  with  a  number  of 

eigenvalues  Ap ,  where  p  =  1 ,2,3 .  For  a  given  frequency,  the  wave  equation  in  the 

waveguide  determines  the  wavenumber  kp  for  each  value  of  p , 

kp  =  to  tc  -  Ap  . 
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When  the  frequency  is  below  cut-off,  oj/c<Ap,  the  wavenumber  k  is  zero,  or 
imaginary,  and  mode  does  not  propagate;  top  =  c  Ap  is  called  the  cut-off  frequency  of  the 
waveguide. 

All  waveguide  modes  have  the  form  £_)_«  e±ikpz-im  The  packwarcj  propagating 

__  ^  ^ 

wave,  E^e  p  ,  results  in  an  FEL  interaction  at  rather  long  wavelengths.  It  can 
sometimes  be  of  interest,  but  most  applications  of  the  relativistic  FEL  interaction  seek 
shorter  wavelengths.  Therefore,  the  backward  propagating  wave  is  not  discussed  here, 
and  we  concentrate  on  waveguide  modes  with  form  E  ^  e+lkpZ~liSit .  When  the 
polarization  of  the  waveguide  mode  is  chosen  to  match  the  electron  motion  in  the 
periodic  undulator  field,  we  have  the  best  coupling.  The  form  of  the  transverse  electron 
motion  is  0^°=  eik°z  where  A0=2rD'/fe0  is  the  undulator  wavelength.  The  fourth 
component  of  the  Lorentz  force  equation  governs  the  electron  energy  evolution  and 
bunching.  It  has  the  form  yocfy^oc  e‘(ko+kp)z -not ,  ^  naturally  defines  the  electron 
phase  £  =  (k0+kp)z-wt  where  k0,  kp,  and  to  are  fixed  by  the  FEL  design,  and  £(r )  follows 
the  evolution  of  z(r).  The  corresponding  the  phase  velocity  is  ~  £.  Generally,  the 
electron  phase  and  phase  velocity  evolve  over  the  whole  undulator  length  L  so  that  it  is 
natural  to  relate  the  time  t  to  the  evolution  time  L/c(3z  = Lie  in  the  relativistic  FEL  with 
electron  beam  z  velocity  (3Z .  The  natural  definition  for  the  electron  phase  velocity  is  then 

‘  r  g  *v 

v  =  L[(k0+kp)Pz-(d/c]=  L  k0+  -A„2  Pz-—  • 

LI  L  c  J  J  c . 

The  electron  phase  velocity  v(r)  depends  on  kQ,  kp,  oo,  and  L,  which  are  fixed  in  an  FEL 
design,  and  follows  the  evolution  of  the  electron  z  velocity  (3z(r). 

The  electron  motion  in  the  periodic  undulator  and  interaction  with  the  waveguide 
mode  are  resonant  when  v  =  0.  This  occurs  at  the  resonant  frequency  to*  =c(k0+kp)$: 
where  kp  =(w*2/c2 -  Ap)V2.  This  can  solved  for  co*  in  terms  of  k0 ,  kp,  and  (3Z.  The 
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eigenvalue  Ap  depends  on  the  waveguide  dimension  and  shape.  For  the  rectangular 
waveguide  with  sides  a  and  b  we  have 


r  2  2 

m  n 

A__  =  K  — —  +  — 
2  .2 
la  6 


where  m,n  =0,1, 2, 3,....,  but  not  both  m  =  n  =  0.  For  the  circular  waveguide  with  radius 


A„„  =x„„/a  , 

nm  run  1 


where  is  the  root  of  J/ix^)  =  0.  Looking  at  the  first  few  roots,  we  see  that  roughly, 
x*m  =nm-  For  both  waveguides,  the  cut-off  frequency  is  cop  =  c  Ap .  The  NRL  waveguide 
is  circular  with  a  =  1.5cm  so  that  the  lowest  values  are  roughly  Ap  ~nm/a  and  give  the 
cut-off  frequency  (op  =  cnm/a  =  10GHz.  The  NRL  FEL  operates  above  cut-off  at 
co  =  35GHz.  Far  above  cut-off,  the  expression  for  the  electron  phase  velocity  can  be 
expanded  in  Ap .  This  give  the  lowest  order  shift  in  phase  velocity  and  clarifies  the  effect 
of  the  waveguide  and  the  FEL  performance. 

Expanding  in  Ap  gives  the  form  of  the  phase  velocity  v  =  v(0> +  Avp  +  ....  where  v(0> 
is  the  old  phase  velocity  definition  without  waveguide  corrections, 
v(0)  =  L[(fc0+co/c)P2-(o/c],  and  Avp  is  the  first  order  waveguide  correction  in  Ap.  The  old, 
unperturbed  resonance  condition,  v(O)  =  0  gives  the  resonant  frequency 
co*  =ck0/(\-$z)~  2Y2cfc0/(1+K2).  A  relativistic  FEL  that  is  near  resonance  co-co*  for 
maximum  coupling,  has  a  shift  in  phase  velocities  caused  by  the  waveguide, 


Avp  ~  - 


A^o2A2(1+*2) 


where  N  =L/X0  is  the  number  of  undulator  periods.  The  shift  in  resonance  is  negative, 
and  diminishes  when  y  -» «>.  in  a  rectangular  waveguide 
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kN\%C\+K?-) 

8y2 


As  an  example,  consider  the  LLNL  ELF  experiment  where  y  =  7,  Xo  =  10cm,  N  =30, 
/;  =2.4,  and  waveguide  dimension  is  b=  3cm.  The  waveguide  shift  is  then 
Av0ELF  =  —  6 7c n2,  and  AvE,LF  = -6k  for  the  lowest  order  mode.  This  is  a  relatively  large 
shift  should  be  detectable.  The  gain  spectrum  bandwidth  for  ELF  is  4/1/6  =  4rc  is 
comparable. 

The  NRL  FEL  uses  a  circular  waveguide;  the  waveguide  shift  in  a  circular 
waveguide  is 


8rcy2 


l  a  J 


For  the  NRL  parameters  y=2.4,  X0  =  4cm,  N~  22,  and  a  =  1.5cm,  the  shift  is 
A^l  ~-Km2/2,  and  a^rl=-tc/2  in  the  lowest  order  mode.  The  gain  spectrum 
bandwidth  is  4 j  ~2 n  so  that  several  modes  could  be  near  resonance  with  good 
coupling. 


Gain  Degradation  and  Electron  Beam  Quality  (Summary  of  Appendix  A) 

Maintaining  the  coherence  of  the  electron  bunches  over  a  significant  interaction 
length  imposes  important  restrictions  on  the  electron-beam  quality.  An  energy  or 
angular  spread  (due  to  emittance)  in  the  NRL  experiment  would  contribute  a  random 
component  to  the  electron  motion  that  decreases  the  coherent  bunching  in  time. 
Accelerators  present  a  design  trade  off  between  high-current  and  high-beam  quality  that 
makes  it  essential  to  accurately  evaluate  the  effects  of  beam  quality  in  present  and 
future  experiments.  The  theoretical  approach  presented  in  the  publication  of  Appendix 
A  uses  a  convenient,  yet  powerful,  method  of  including  an  arbitrary  electron  distribution 
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function  in  a  self-consistent  integral  equation  for  the  complex  optical  field.  The  FEL  is 
described  by  solving  the  Lorentz-Maxwell  equations  self-consistently  in  weak  optical 
fields.  The  field  evolution  is  determined  by  an  integral  equation  that  allows  the  inclusion 
of  an  arbitrary  electron  distribution  function  in  a  simple  way.  FEL  gain  and  the  effects  of 
beam  quality  can  then  be  calculated  analytically  or  integrated  on  a  small  computer.  In 
strong  optical  fields,  the  effects  of  beam  quality  are  diminished  because  large  optical 
buckets  induce  and  even  larger  spread  than  in  the  initial  beam. 

Computer  simulations  are  not  as  useful  in  solving  this  kind  of  FEL  problem  because 
they  use  a  prohibitively  large  number  of  sample  particles,  and  introduce  a  large  amount 
of  numerical  noise  when  distributed  over  a  large  volume  of  phase  space. 

In  Appendix  A,  contour  maps  are  used  to  show  the  gain  degradation  due  to  an 
electron-beam  energy  spread  and  an  electron-beam  angular  spread.  Figures  3  and  4 
are  specifically  chosen  to  represent  the  NRL  experiment.  They  show  that  the  amount  of 
beam  spread  expected  in  the  NRL  experiment  can  diminish  gain.  Furthermore,  the  two 
distributions  give  quite  different  results  for  the  same  amount  of  spread  emphasizing  the 
importance  of  the  detailed  shape  of  the  distribution  function.  When  there  is  high  gain,  as 
in  the  NRL  experiment,  the  exponential  growth-rate  becomes  less  susceptible  to 
degradation  from  the  electron-beam  quality.  Other  types  of  distribution  functions  can  be 
added  to  the  theoretical  technique  developed.  If  solved  numerically,  even  experimental 
distribution  functions  peculiar  to  a  given  accelerator  or  transport  system  can  be  added. 

This  research  resulted  in  the  publication  included  in  Appendix  A.  The  reference  is 
W.  B.  Colson,  J.  C.  Gallardo,  and  P.  M.  Bosco,  "Free-Electron  Laser  Gain  Degradation 
and  Electron  Beam  Quality",  Physical  Review  A34,  4875  (1986). 
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Electron  Trapping  in  FEL  Oscillators  and  Amplifiers  (Summary  of  Appendix  B) 

In  high-power  FELs,  the  electrons  can  become  trapped  in  deep  potential  wells 
formed  by  the  combined  optical  and  undulator  field  forces.  The  trapped  current 
oscillates  at  the  synchrotron  frequency,  and  can  drive  the  optical  wave  at  sideband 
frequencies  around  the  fundamental.  This  "trapped-particle  instability"  can  occur  in  both 
the  oscillator  and  amplifier  configurations.  This  is  an  undesirable  effect  for  many 
applications,  and  would  certainly  be  a  serious  problem  for  the  two-stage  FEL. 

The  trapped-particle  instability  starts  when  electrons  become  trapped  in  deep 
potential  wells  formed  by  the  combined  optical  and  undulator  fields.  Electrons  near  the 
bottom  of  the  well  oscillate  in  harmonic  orbits  at  the  synchrotron  frequency  causing  the 
carrier  wave  to  develop  sidebands.  The  instability  has  been  clearly  observed  in  the 
high-power  FEL  oscillator  at  LANL,  and  possibly  in  the  TRW/Stanford  oscillator.  Here, 
the  multimode  sideband  theory  is  reviewed  with  examples  from  oscillators,  and 
amplifiers.  Common  features  and  differences  are  discussed  along  with  there  relevance 
to  the  NRL  experiment. 

Increasing  the  current  density  j  or  the  resonator  Q  increases  the  steady-state 
power,  the  synchrotron  frequency,  and  the  sideband  gain.  Over  many  hundreds  of 
passes,  the  stored  optical  wave  "sees"  many  synchrotron  oscillations,  so  that  any 
sideband  gain  above  threshold  gives  large  growth  from  a  small  amount  of  noise.  The 
resulting  steady-state  features  are  therefore  not  affected  by  the  details  of  the  noise 
source  employed.  The  addition  of  sideband  power  is  cumulative,  since  the  presence  of 
a  strong  sideband  again  increases  the  steady-state  power.  This  gives  a  broad,  possibly 
chaotic,  optical  power  spectrum,  and  a  broad  electron  distribution.  If  taper  is  introduce 
into  the  undulator  design,  the  synchrotron  frequency  is  only  slightly  modified,  and  the 
sideband  gain  is  reduced. 

If  the  electron  pulse  is  not  short,  it  is  prudent  to  simulate  the  FEL  by  sampling  a 
smaller  window  with  periodic  boundary  conditions.  In  the  long-pulse  limit,  as  in  the  NRL 


experiment,  the  wrapped-window  approach  is  the  only  viable  solution  to  the  FEL 
multimode  simulation. 

In  the  high-current  FEL  amplifier,  y»1,  large  optical  power  can  be  produced  in  a 
single  pass  through  the  undulator.  The  FEL  growth  rates  are  so  large  that  electrons 
become  trapped  early  in  the  undulator  and  begin  executing  synchrotron  oscillations. 
Even  the  limited  number  of  synchrotron  cycles  can  result  in  significant  sideband  gain 
owing  to  the  large  current  density  j.  The  trapped-particle  instability  in  FEL  amplifiers 
differs  in  many  ways  from  the  oscillator  case.  Sources  of  noise  are  much  more 
important,  and  there  is  no  resonator  Q  to  consider.  When  the  sideband  instability  does 
occur,  the  FEL  does  not  reach  steady-state  operation  even  in  strong  optical  fields.  The 
power  continues  to  increase  and  so  does  the  synchrotron  frequency. 

This  research  resulted  in  the  publication  included  in  Appendix  B.  The  reference  is 
W.B.  Colson,  "The  Effect  of  Electron  Trapping  in  Free-Electron  Laser  Oscillators  and 
Amplifiers",  Proceedings  of  the  1985  International  Conf.  on  LASERS,  Las  Vegas  CA 
(1986). 

NRL  FEL  Simulations 

Some  simulations  are  shown  below  that  correspond  to  the  parameters  of  the  NRL 
FEL  experiment.  In  the  next  figures,  the  current  density  j  =100  is  used  with  Q  =2  and 
<3=3.  Details  of  the  figures  are  explained  in  Appendix  B.  The  evolution  of  the  optical 
field,  optical  spectrum,  electron  spectrum,  net  gain  G,  and  power  are  shown.  Each 
simulations  starts  from  shot  noise  and  runs  for  n  =  200  passes.  Both  cases  show  that 
the  sideband  instability  has  increased  the  optical  spectrum  P{v)  beyond  the  gain 
spectrum  bandwidth  in  weak  fields  G(v).  When  the  Q  is  increased,  the  power 
increases,  but  the  spectral  quality  decreases.  The  field  amplitude  |  a  (z  ,n )  |  shown  in 
the  upper  left  corner  would  provide  the  electromagnetic  undulator  in  the  two-stage  NRL 
FEL.  Because  of  the  sideband  instability,  it  would  be  a  poor  quality  undulator. 
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On  the  previous  page  the  figures  show  simulation  results  with  the  current  density 
j  =300  using  Q  =  3  and  <2=4.  The  simulations  are  started  from  shot  noise  and  run  for 
200  passes  as  in  the  first  two  figures.  The  increased  current  increases  the  final  field 
strength  dramatically  from  about  |  a  \  ~  250  to  |  a  |  -  900.  The  increased  power  spreads 
the  electrons  over  a  wider  range  of  phase  velocities,  but  also  degrades  the  optical 
spectrum  />(v).  The  optical  field  and  spectrum  are  now  chaotic  and  broad  band. 
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ABSTRACT 


The  free  electron  laser  can  be  described  by  solving  the  Lorentz-Maxwell  equations 
self-consistently  in  weak  optical  fields.  The  field  evolution  is  determined  by  an  integral 
equation  that  allows  the  inclusion  of  an  arbitrary  electron  distribution  function  in  a 
simple  way.  Contour  maps  are  used  to  show  the  gain  degradation  due  to  an  electron 
beam  energy  spread  and  an  electron  beam  angular  spread.  In  the  limit  of  low  gain,  the 
gain  spectrum  is  related  to  the  spontaneous  emission  line-shape  through  successively 
higher  derivatives.  In  the  limit  of  high  gain,  it  is  shown  that  the  growth  rate  becomes 
less  susceptible  to  degradation  from  the  electron  beam  quality. 

I.  Introduction 

In  a  tree-electron  laser  (FEL),  a  relativistic  electron  beam  amplifies  a  co-propagating, 
coherent  optical  wave  traveling  through  a  periodic  undulator  magnetic  field  [1],  In  the  oscillator 
configuration,  coherent  electron  bunching  develops  on  each  pass  while  resonator  mirrors  allow 
the  stored  optical  power  to  grow  over  many  passes.  In  the  amplifier  configuration,  coherent 
electron  bunches  develop  rapidly  in  the  first  part  of  the  undulator  followed  by  rapid  growth  of  the 
optical  field.  Maintaining  the  coherence  of  the  electron  bunches  over  a  significant  interaction 
length  imposes  important  restrictions  on  the  electron  beam  quality.  An  energy  or  angular  spread 
(due  to  emittance)  contributes  a  random  component  to  the  electron  motion  that  decreases  the 
coherent  bunching  in  time. 
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Some  of  the  earliesl  FEL  experiments  used  electron  beams  that  were  essentially 
monoenergetic  [2-5],  but  practically  all  subsequent  experiments  have  made  use  of  higher  current 
sources  with  significant  energy  spread  or  emittance.  Many  accelerators  present  a  design  trade-off 
between  high  current  and  high  beam  quality.  This  makes  it  essential  to  accurately  evaluate  the 
effects  of  beam  quality  in  present  and  future  experiments.  It  is  particularly  important  for  FELs 
designed  to  operate  at  XUV  or  X-ray  wavelengths  [6].  Several  theoretical  models  involving 
simulations  and  plasma  dispersion  relations  have  discussed  the  detrimental  effects  of  electron 
beam  quality  in  the  FEL  interaction  [7-24].  The  theory  presented  here  uses  a  convenient,  yet 
powerful,  method  of  including  an  arbitrary  electron  distribution  function  in  a  seff-consistent  integral 
equation  for  the  complex  optical  field.  FEL  gain  and  the  effects  of  beam  quality  can  then  be 
calculated  analytically  or  integrated  on  a  small  computer. 

Since  the  basic  equations  solved  here  are  the  same  as  in  computer  simulations  or  the 
plasma  dispersion  methods,  specific  physical  results  have  been  shown  to  agree  with  those 
methods  when  a  direct  comparison  is  possible.  The  computer  simulations  have  proved  to  be  a 
useful  method  of  understanding  many  aspects  of  the  FEL  interaction,  but  one  of  the  most  difficult 
effects  to  accurately  characterize  is  that  of  electron  beam  quality.  Even  a  prohibitively  large 
number  of  sample  particles  is  far  short  of  the  number  in  a  real  experiment,  and  yet  introduces  a 
large  amount  of  numerical  noise  when  distributed  over  a  large  volume  of  phase-space.  To 
reproduce  some  of  the  results  shown  later  in  this  paper,  we  found  the  simulation  method  to  be 
several  hundred  to  a  thousand  times  less  efficient.  While  many  other  FEL  topics  are  most 
efficiently  studied  through  simulations,  the  detrimental  effects  of  beam  quality  are  probably  better 
handled  through  a  combination  of  analytic  and  numerical  techniques.  The  stability  analysis  used 
to  obtain  plasma  dispersion  relations  usually  calculates  the  reduced  FEL  growth  rates  due  to  poor 
beam  quality.  This  method  can  lead  to  analytical  expressions,  but  depends  upon  specific  models 
for  the  electron  beam  distribution,  and  does  not  easily  describe  more  complicated  transient 
behavior  where  the  FEL  growth  rate  is  not  constant;  the  FEL  is  often  designed  to  operate  in  this 
regime.  In  addition,  the  exact  formulation  presented  here  works  smoothly  between  different 
regimes  of  operation  like  high  and  low  gain.  The  only  requirement  is  weak  optical  fields. 
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II.  Basic  Theory 

We  solve  the  electron  Lorentz  and  optical  wave  equations  self-consistently  with  the 
assumption  of  weak  optical  fields.  The  effects  of  beam  quality  are  typically  less  important  when 
the  optical  field  strength  is  large  near  saturation,  and  the  issue  of  beam  quality  is  most  important 
in  weak  fields  where  the  accurate  evaluation  of  gain  can  determine  whether  the  FEL  is  above  or 
below  threshold. 

The  electrons  travel  through  a  periodic  undulator  with  the  field  on  the  z  -axis  described  by 
i?=B  [cos(Jfcoz).sin(jfc0z),0]  where  B  is  the  peak  magnetic  field  amplitude.  The  undulator  field 
extends  over  a  length  L  =  NX0  with  a  number  of  periods  N ,  and  wavelength  Xo  =  2jt /k0.  The 
electron  velocity  in  a  perfect  helical  orbit  is  cj?  =  c  [-(K/y)  cos(jfc0z),-  (Kly)  sin(*0z),  p0]  where 
K  =  eB\0i2nmc2,  e  is  the  electron's  charge  magnitude,  m  is  the  electron's  mass,  c  is  the  speed  of 
light,  p0  =  ( 1-(1  +K  2)/y2  )v2,  and  ymc2  is  the  electron's  energy.  Imperfect  injection  due  to  poor 
beam  quality  is  more  meaningfully  introduced  after  some  further  theoretical  development.  A 
typical  undulator  uses  B  =  2kG  and  Xq  =  5 cm,  so  that  K  =  1.  Since  the  electrons  are  relativistic 
(y  »1),  the  transverse  excursions  are  small  compared  to  V 

The  optical  vector  potential  with  the  polarization  that  best  couples  to  the  above  trajectory  is 
A  -  fc-1 1  £  |  [  simy,  cosvy,  0  ]  where  \y  =  *z-GM+q,  and  X  =  2nlk  =  2nc  /to  is  the  optical  carrier 
wavelength.  The  complex  electric  field  envelope,  £(z,r)=|£(z,r)|e‘*<z'f)  ,  is  taken  to  vary  slowly 
in  z  and  t ,  so  that  terms  containing  two  derivatives  in  the  wave  equation  are  small  compared  to 
terms  with  single  derivative  [25].  No  transverse  (x.y)-dependence  is  included  so  that  diffraction  is 
taken  to  be  a  small  over  the  interaction  length  L ,  and  the  electron  beam  remains  aligned  near  the 
center  of  the  optical  mode.  The  transverse  motion  above,  proportional  to  {K/y),  defines  the 
transverse  current  for  each  electron  in  the  beam.  If  the  current  density  is  uniform  over  a  sufficient 
length,  each  point  z  +  ct  in  the  optical  field  envelope  evolves  according  *o  the  slowly-varying  wave 
equation  [25] 


where  a  =  4NneKLE/y2mc2  is  the  dimensionless  optical  field  strength,  x  =  ct/L  is  the 
dimensionless  time  (0<x<  1),  j  =8N(neKL)2p /y2mc2  is  the  dimensionless  current  density,  p  is 
the  actual  electron  particle  density,  £  =  ( k+k0)z -  au  is  the  electron  phase  in  the  combined  optical 
and  undulator  fields,  and  <...>  represents  a  normalized  average  over  all  electrons  in  the  beam 


driving  a(x).  The  electrons  are  labeled  by  their  initial  phase-space  coordinates;  the  initial  phase  is 
=  £(0)  ,  and  the  initial  phase  velocity  is  \i=dC,(0)/dx  =  L  [(k  +k0)%- k] .  There  are  a  large 
number  of  electrons  spread  randomly  over  each  optical  wavelength  (-107 ),  so  that  the  C,  can  be 
accurately  taken  to  be  uniformly  spread  along  each  section  of  the  electron  beam  one  wavelength 
of  light  long.  It  can  be  easily  seen  in  (1)  that  bunching  the  electrons  near  the  relative  phase 
£  +  <(>  =  n  drives  the  optical  wave  amplitude  producing  gain,  while  bunching  near  £  +  q  =  rc/2  drives 
the  optical  phase  <)>  without  gain.  Bunching  electrons  near  £  +  q  =  0  results  in  negative  gain,  or 
absorption.  The  dimensionless  electron  phase  velocity  v,  has  an  initial  spread  associated  with  the 
beam  quality. 

The  electron  motion  in  the  presence  of  the  optical  wave  is  described  by  the  Lorentz  force 
equation;  dyldt  =  -(eimc)#-I?.  In  the  FEL,  it  is  important  to  distinguish  between  collective 
Coulomb  forces  and  collective  high-gain  effects  [26].  Most  FELs  do  not  use  current  densities 
large  enough  for  Coulomb  forces  to  be  a  significant  effect  for  the  relativistic  electrons;  yet,  high 
gain  is  possible  and  will  be  included.  Using  the  definitions  and  assumptions  above,  the  Lorentz 
force  takes  on  the  form  of  the  pendulum  equation  [27], 


dzt  dv 

— -  =  —  =  |  a  |  COS  (C+4>) 
dxz  dx 


(2) 


The  combined  equations  (1)  and  (2)  are  valid  in  weak  or  strong  optical  fields,  for  large  or  small 
gain,  and  for  an  arbitrary  electron  distribution.  Strong  fields  near  saturation  mean  that  [a  |  »  rc, 
and  weak  fields  occur  when  \a  \  «n.  High  gain  is  achieved  when  j  »  1  and  low  gain  occurs 
when  j  <  1  [28].  Useful  FEL  configurations  display  a  wide  range  of  current  densities.  The 
electron  beam  area  is  typically  between  1mm  and  5mm,  but  the  current  ranges  from  1A  up  to 
lOkA.  Undulator  lengths  L  now  range  from  1m  to  5m,  but  will  soon  be  made  to  L  =20m  and 
beyond.  With  electron  energies  in  the  range  lOMeV  to  iGeV,  the  corresponding  values  of  j  are 
from  unity  to  more  than  5x1 04  [1].  Both  the  high-gain,  single-pass  and  the  low-gain,  oscillator 
configurations  have  important  applications. 

Equations  (1)  and  (2)  were  originally  derived  [25]  for  the  more  general  case  where  the 
electron  energy  can  change  significantly  during  a  single  pass;  in  this  case,  an  additional  factor 
r|  =  (1  -  v/2kN)  alters  the  wave  and  electron  equations  so  that  d  =  -j  <  e~'t'  >  and, 

C  =  v  =  | a  |ti  2cos  (£  +  <+>)  with  (')  =d(  )!dx  .  The  following  work,  however,  will  be  confined  to  weak 
optical  fields  where  r\  ~  1.  An  extension  to  higher  harmonics  and  linearly  polarized  undulators  is 
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also  possible  without  any  change  in  form  of  (1)  and  (2),  so  that  the  general  conclusions  and 
methods  of  this  paper  are  immediately  applicable  to  a  wide  range  of  FEL  designs. 

We  now  proceed  to  solve  (1)  and  (2)  in  weak  fields,  \a  \ «  k,  to  obtain  an  integral  equation 
for  a(x)  incorporating  an  arbitrary  electron  distribution  function.  The  electron  phase  can  be 
expressed  as  C  =  ^  +  v;t  +  £ (1)  where  £ !1)  is  the  first-order  perturbation  in  a .  Expanding  (1)  and 
(2)  we  have 

t 

a(x)  =  a0  +  ijjQds<  exp[-i(  £  +  v;i  )  ]C(1)(^)  >  -  (3) 

-J  1  q 

C(1>(s)  =  —jodqjodu[a(u)  exp(  i(^ +  vtu  ))  +  a*  (u)  exp(  -i(  £,•  +  v;u  ) )  ]  , 
where  the  initial  optical  field  is  a(0)  =  |a(0)j  =a0  and  <f>(0)  =  0.  We  have  made  use  of 

2k 

<exp(-«C,  )>=  Jq  ^^exp( -r'C  )/2k  =  0  ,  since  the  initial  electrons  are  spread  uniformly  in 
phase.  The  reference  to  the  individual  electron  phases  C(1>  can  be  explicitly  removed  by 
combining  the  equations  in  (3).  Then,  we  have  an  integral  equation  governing  the  evolution  of  the 
optical  field  a  (t): 

ij  ^  ^  ? 

a(x)  =  aQ  +  -^jodsjQdqjodu  <exp(-ivi(s-u))>a(u)  ,  ,  (4) 

where  <...>  is  now  an  average  over  the  initial  electron  velocity  distribution,  and  all  reference  to  the 
electron  phases  has  been  removed.  Since  (4)  is  an  iterated  triple  integral,  it  may  be  rewritten  [29] 
as  a  double  integral, 

ij  x  1 

a(*)~aQ  +  —j^ds^dq  <  exp  (-iVi(s-q))  >  (s-q)  a{q)  .  (5) 

A  normalized  electron  distribution  function  /  (v()  can  be  used  to  evaluate  the  remaining  average: 

*•  M 

<  > 3LdV; /(vi)(-) with L^/(v.) = i- 

HI.  Simple  Electron  Distributions 

We  begin  by  considering  two  simple  examples  with  perfect  beam  quality.  In  the  first,  we 
start  the  FEL  on  resonance  where  the  electron-optical  wave  coupling  is  largest,  /(v,)  =  5(vI). 
The  optical  wave  is  most  simply  determined  from  (4). 


(6) 
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ij  *  1  tq 

a(x)  =  a 0+  —jQdsjodqj^du  a(u) 

The  integral  equation  (6)  can  also  be  written  in  a  differential  form  by  taking  successive  derivatives, 

3 

a  (x)  =  ija{x)/2.  The  complete  solution  uses  the  form  a (x)  =  £a„exp(  a** )  where  the  aH  are  the 

«=1 

three  complex  roots  of  the  cubic  equation  a3  -  iji2  =  0,  and  the  coefficients  an  are  determined  by 
the  initial  conditions  a (0)  =  a0 ,  anda(O)  =  a(0)  =  0  [1J.  The  solution  lor a(x)  is 

a (x)  =  —  I  exp  ( ( j/2  )1/3  M 3 )  x/2  )  +  exp  ( ( j/2  )1/3  (i-V 3  )xl2)  +  exp  (~i(jl2 )1/3 x )  .  (7) 

3  L 

If  the  current  density  is  small  j  ->0,  orx  « 1  ,  we  have  the  trivial  result  a(x )  =  a0(t+ijx3lt2+... ) 
for  an  FEL  starting  on  resonance.  There  is  no  change  in  the  optical  amplitude  J  a  |  =  a0+...  to 
lowest  order,  and  therefore  no  gain.  The  optical  phase  q(x)  increases  slowly  in  proportion  to  t3. 
The  FEL  gain  is  defined  as  G  (x)  =  (  |  a  (x)  | 2  -  a  l  )/a  |  .and 

G  (x)  =  ^  [  2  cosh  ( ( j/2  )1/3  <3  X )  +  4  COS  ( ( ,72  )1/3  3  x/2 )  COSh  ((j/  2  ) 1/3  ^3  x/2  )  -  6  ]  .  (8) 

In  the  high  current  limit,  ;  »1  on  resonance,  the  expressions  simplify  because  one  fastest- 
growing  root  dominates  and  describes  exponential  growth  in  x.  As  seen  from  (7)  there  is  little 
change  in  the  field  dunng  the  bunching  iime,  x  <  xB  -  ( 2 Ij  )  ,  that  precedes  exponential  growth. 
During  this  time,  the  electrons  move  from  their  initially  uniform  phase  distribution  to  bunch  near 
the  phase  £  +  <j>  =  rc/2.  As  soon  as  bunching  forms,  the  high  current  immediately  causes 
exponential  field  growth  and  high  gain.  Then, 

a(x)  =  —  exp[(;/2)1/3V3T/2]  ,  and  G(x)  =  —  exp[  (  j/2  )1/3^3x]  .  (9) 

3  9 

A  second  simple  example  is  a  high  quality  electron  beam  starting  off  resonance  at  v0.  This 
is  characterized  by  /  (v,)  =  5(  v,  -  v0 ).  The  optical  field  is  then  determined  by 

a(x)  =  a0  +  ^\Qds^dqjQdue  ‘V,U'“ 1  a{u)  .  (10) 

For  low  current,  j<  1,  the  optical  field  evolution  away  from  a0  is  small  so  that  a(u)  =aQc an  can  be 
extracted  from  the  integrand  in  (10).  The  resulting  integrals  are  easily  solved  to  obtained  the 
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usual  low  current  gain  and  phase  shift  formulas  [27], 

2  -  2cos(v0t)  -  v0t  sin(v0-c)  2  sin(v0x)  -  v0t  (  1  +  cos(v0x) ) 

G  (x)  =  ;' - - -  ,  and  4»(t)  =  j - - -  (11) 

Vo  2v03 

The  maximum  final  gain  is  G  =0.135;  at  v0  =  2.6  and  x  =  1 ,  while  the  range  of  modes  with 
significant  gain  is  8v0  =  1  about  the  peak. 

In  order  to  obtain  the  general  solution,  use  the  substitution  b  =a  <?‘V°T  in  (10).  Successive 
derivatives  then  lead  to  the  differential  form  of  (10),  b  -i\Jb  -ijbl2  =  0.  Solutions  of  the  form 

3 

b  =  £b„exp(  aBx )  have  roots  a,  that  satisfy  the  cubic  equation  a3-iv0a3-t;72  =  0.  In  the  limit 

of  high  current ;'  »  1 ,  the  exponential  gain  coefficient  is  reduced  by  the  factor  -  v|  /  3^3  ( ji 2  )1/3 
so  that  the  gain  spectrum  is  centered  about  v0  =  0  with  a  characteristic  range  8v0  =  4.22; 1/6.  In  the 
high  current  case,  the  range  of  modes  with  significant  gain  increases  slowly  as  j  increases.  We 
go  on  now  to  look  at  more  interesting  FEL  distributions  describing  less  than  perfect  electron 
beams. 

IV.  More  General  Electron  Distributions 

New  cases  of  interest  involve  more  complicated  distributions  / (v,)  describing  the  initial 
electron  beam  in  the  integral  equation  (5).  Two  electrons  starting  at  the  same  phase  C,  at  the 
beginning  of  the  undulator  (x  =  0),  but  with  different  z-velocities,  c  p0  and  c  (P0  +  ^Po).  will  drift  apart 
as  they  travel  through  the  undulator.  The  amount  of  drift  is  not  easily  predicted  without  solving  the 
full  problem,  because  electrons  can  influence  each  other  through  the  self-consistently  evolving 
optical  wave.  In  this  sense,  the  effect  of  FEL  beam  quality  is  collective.  However,  the  times  for 
the  two  electrons  to  traverse  the  undulator  are  nearly  identical  since  they  are  relativistic,  Li  p0c  = 
L/(p0  +  Ap0)c  ~Uc.  An  estimate  of  their  separation  at  the  end  of  the  undulator  (t=1)  is 
A z  =ApoL,  and  their  approximate  phase  difference  is  A (,=  (k  +  k0)Az  =kAz  =kL Ap0  •  If  the 
velocity  difference  Ap0  is  due  to  an  initial  energy  difference  A ymc2 ,  we  have  Ap0  =  (1+/C  2)Ayi  y3 , 
and  an  approximate  final  phase  separation  A  £  =  4n w  Ayly. 

Any  random  phase  difference  A£  -  rt,  or  larger,  between  electrons  in  the  beam  is  important 
to  the  FEL  operation,  because  bunching  on  the  optical  wavelength  scale  is  diminished 
significantly.  At  the  end  of  the  undufafoi,  the  final  phase  difference  is  roughly  estimated  by 
AC  =  AVj  for  each  electron.  From  the  definition  of  the  electron  phase  velocity  v,  ,  we  see  that  a 


small  change  in  the  initial  electron  energy  Ay  me 2  corresponds  to  a  change  in  the  initial  phase 
velocity,  Av,  =  4?&VAy/yfory»  1.  A  distribution  of  initial  electron  energies  from  an  accelerator  or 
storage  ring  is  often  accurately  represented  by  the  normal  distribution  function  so  that  we  can  take 


/(V;)  =  - 


exp(-  (v,  -  v0) 2 /2a2) 


where  a  is  the  standard  deviation  of  v,  away  from  the  peak  phase  velocity  v0.  If  Ay  me 2  is  taken  to 
be  the  standard  deviation  of  the  electron  energy  away  from  y me2,  then  cr  =  47tNAy/y  .  Two 
electrons  starting  at  the  same  phase  ^  ,  but  with  an  energy  difference  Ay=  y  /  4N  will  drift  apart  by 
roughly  half  of  one  optical  wavelength  at  the  end  of  the  unduiator.  A  random  spread  of  width  a  =  k 
causes  a  random  phase  spread  of  approximately  a£  =  n  at  the  end  of  the  unduiator  and  impairs 
bunching.  Inserting  (12)  into  (5)  gives 

a(x)=a0+^fodsfodq  (s-q)a(q)  .  (13) 

The  gaussian  factor  in  the  integrand  decreases  the  coupling  current  j  as  t  increases,  and 
describes  the  degradation  of  bunching  due  to  the  spread  in  electron  phase  velocities.  The 
complicated  self-consistent  evolution  of  the  electron  beam  distribution  and  the  optical  field  are 
described  exactly  in  (13),  but  before  evaluation,  we  can  generalize  its  form  further. 

An  angular  spread  is  also  possible  due  to  the  finite  emittance  of  an  electron  beam.  An 
electron  of  energy  ymc2  entering  the  unduiator  with  a  small  injection  angle  9;  has  a  reduced 
velocity,  po-»pocos0;  =p0(1-9,2/2)  .  The  resulting  z -velocity  change  is  Ap0  =  -92/2  , 
reducing  the  initial  phase  velocity  by  Av,  =  -  27t/Vy292/(1+A: 2 )  .  A  gaussian  distribution  of  angles 
about  the  z  -axis  with  standard  deviation  A9;  gives  the  exponential  distribution  function 


/(v;)  =  - 


exp  (  -  ( v0  -  V,- )  /  o9 ; 


for  v,<v0  ,  and  /(v;)  =  0  for  v,>v0  ,  (14) 


where  <j0=  4itA/y2A9,  Z/(1+K  2),  and  v0  is  the  phase  velocity  for  electrons  entering  on-axis.  The 
distribution  function  (14)  is  sharply  peaked  at  v,  =  v0  where  electrons  enter  on-axis,  and  decays 
exponentially  for  v<  <  v0  because  the  injection  at  any  angle  9,  can  only  decrease  the  electron's 
longitudinal  velocity  and  its  phase  velocity  v,.  If  each  element  of  the  energy  distribution  (12)  is 
given  an  angular  spread  according  to  (14),  then  the  resulting  integral  equation  for  the  optical  field 
becomes 


(  s-q)  a(q) 


(15) 
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a(x )  =a0  + 


1  -  i  Gg(  s  -q ) 


-<v0(  s-q) 

e 


The  transverse  motion  of  electrons  injected  at  an  angle  is  sometimes  confined  by  either  the 
natural  off-axis  undulator  fields  or  external  focusing  elements.  The  focusing  forces  result  in 
transverse  betatron  oscillations  about  the  undulator  axis.  When  the  electron  beam  is  injected  to 
match  the  natural  focusing  properties  of  an  undulator,  the  number  of  betatron  oscillations  along 
the  undulator  is  N^  =  NK/^2y.  In  the  limit  of  large  y  and/or  small  K ,  the  angular  spread  of 
electrons  can  be  important,  cr9  >  1 ,  while  the  transverse  focusing  can  be  made  negligible,  N p  <  1  . 
In  this  limit,  the  integral  equation  (15)  applies. 

The  complex  optical  field  a(x)  now  depends  on  an  input  energy  spread  characterized  by  o  , 
and  an  input  angular  spread  characterized  by  cr9  .  Other  types  of  distribution  functions  can  be 
added  in  a  similar  way.  If  (5)  is  solved  numerically,  even  experimental  distribution  functions 
peculiar  to  a  given  accelerator  or  transport  system  can  be  added.  The  general  result  (5)  and  the 
specific  example  (15)  are  important  results  of  this  paper.  They  provide  analytic  expressions 
describing  FEL  performance  with  an  arbitrary  electron  distribution  function. 


V.  Low-Current  FELs 

One  of  the  cases  of  general  interest  is  the  low  current  FEL  oscillator.  Radiation  energy  is 
stored  in  an  optical  resonator,  and  repeatedly  driven  by  successive  electron  pulses  from  an 
accelerator  like  a  linac  or  storage  ring.  An  important  issue  for  the  oscillator  is  the  detrimental 
effect  of  the  electron  energy  and  angular  distributions  when  the  oscillator  is  starling  from  weak 
optical  fields.  In  the  low-current  case,  we  can  simplify  (15)  by  taking  a(q)  =  a0  in  the  integrand  of 
the  integral  equation,  and  neglecting  higher-order  terms  in  j  .  Without  a(q)  in  the  integrand,  the 
integral  can  be  further  simplified  by  noticing  that  iqe  v’ q  =  -  3Vj  e  v" q .  Then, 


a(x)-a  o 

a0 


j  3  fT  / 
2 


1  -  ('< T9  q 


e 


(16) 


Direct  integration  of  (16)  is  possible,  but  the  result  is  a  complicated  expression  containing  many 
error-functions  [30], 

We  can  alter  the  form  of  (16),  however,  to  obtain  some  important  physical  interpretations. 
The  factors  e"0’2  and  ( 1  —  /  crQ  <7  ) -1  can  be  interpreted  as  power  series  expansions  in  q 


.  Now,  write 


multiplying  the  factor  e-1''1’  ;  then  (-/<? )"  <?  1V’‘7  ->dCte  v,‘l 


a(x)-a0 

a0 


j  exp  (  a  2  32t  /  2  ) 
2  1  +  c8  3Vj 


e 


-*Vj? 


(17) 


The  double  integral  is  simply  v92  ( 1  -  zv0x  -  e~‘v‘z ),  so  that  the  complex  integrations  in  (16)  have 
been  replaced  with  a  power  series  expansion  to  all  orders  in  o  and  c9.  To  first  order  in  a  and  cre, 
an  explicit  expression  for  a(x)  is  easily  obtained  from  (17).  This  is  a  useful  linflSfcg  case  since  a 
low-current  FEL  system  would  not  typically  use  a  low  quality  electron  beam  ( large  a  or  cr9 )  and 
remain  above  threshold.  From  (17),  the  low-current  FEL  gain  at  the  end  of  the  undulator  is 


;  exp(a2a2/2) 

- - a. 

2  1  +  CTg  aVt 


'  sin 2( Vo/2)  ' 
.  (vq/2)2  . 


(18) 


We  recognize  the  factor  in  brackets  [...]  as  the  FEL  spontaneous  emission  line-shape  for  an 
electron  in  a  perfect  trajectory  through  the  undulator.  It  has  been  known  for  some  time  that  the 
gain  is  fundamentally  related  to  the  derivative  of  the  spontaneous  emission  line-shape  [31].  The 
new  feature  presented  in  (18)  is  to  express  how  the  electron  beam  energy  and  angular  spreads 
affect  that  relationship  through  successively  higher  derivatives. 

With  the  physical  interpretation  of  the  line-shape  factor  [...]  ,  we  can  substitute  alternate 
forms.  One  convenient  choice  is  [...]-» exp  (-v02/4tt)  which  approximately  reproduces  the 
correct  features  of  the  simple  FEL  gain  spectrum,  G  =  ( /v0  /4tt  )  exp  ( -vo/4jt ).  The  successive 
derivatives  evaluating  the  effects  of  beam  quality  lead  to  more  compact  expressions,  and  illustrate 
how  (18)  can  be  used  in  practical  situations.  Even  an  experimental  line-shape  could  be  used  in 
(18). 

While  the  analytic  results  presented  have  their  merit,  the  complete  integral  (15)  is  easy  to 
integrate  on  a  small  computer.  The  values  needed  for  the  contour  plots  of  this  paper  were 
evaluated  in  this  way.  Figure  1  shows  a  combined  intensity  and  contour  plot  of  ln(l+G(o,v0)) 
where  the  final  gain  at  the  end  of  the  FEL  undulator  is  G  =(a*(1)a(1)-a|  )/a9  .  In  Fig.  1 
a9  =  0  ,  so  that  gain  degradation  is  only  due  to  an  energy  spread  with  no  angular  spread.  The 
current  density  is  /  =  1 ,  and  gives  low  gain  so  that  ln(1+G)  =  G  .  The  brightest  points  (white)  on 
the  (o,v0)-surface  indicate  peak  gain  G  =0.13/,  while  the  darkest  points  (black)  indicate  maximum 
absorption  G  =-0.13/.  Zero  gain  is  indicated  by  the  intermediate  grey  shown  in  the  scale  at  the 
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top.  Specific  contours  of  constant  gain,  In  (1+G)  =  ±  0.06,  ±0.08,  ±0.1,  and  ±0.12,  are 
superimposed  on  the  intensity  plot.  The  gain  surface  is  approximately  antisymmetric  about  v0  =  0, 
and  in  the  limit  j  -»0,  the  gain  G(c x,v0)  becomes  exactly  antisymmetric.  The  characteristic 
amount  of  spread  required  to  decrease  the  gain  is  seen  from  Fig.  1  and  (18)  to  be  a  *  =  1 .  Note 
also  that  as  the  spread  c  increases,  the  phase  velocity  for  peak  gain  v0*  =  2.6  increases  slightly. 
I  Peak  absorption  occurs  at  -  v0*  and  slightly  decreases  with  increasing  a. 


FEL.  Gain  Surface  1300001 


Fig.  1.  Intensity  and  contour  plot  of  ln(1+G(cx,v0))  with  j  =  1  and  cx9  =  0.  The  weak-field  gain 
degradation  in  this  low-current  FEL  is  due  to  an  electron  beam  energy  spread  with  a  normal 
distribution  function. 

Figure  2  shows  a  combined  intensity  and  contour  plot  of  ln(l+G(a9,v0))  evaluated  by  (15) 
with  j  =  1  and  a  =  0.  The  gain  degradation  here  is  due  to  a  monoenergetic  electron  beam  entering 
the  undulator  with  an  angular  spread  described  by  c0.  The  grey  scale  and  contours  of  gain  are 
the  same  as  in  Fig.  1.  Unlike  Fig.  1,  the  absorption  contours  (white)  have  a  much  different  shape 
than  the  gain  contours  (black).  Since  the  distribution  function /(v,)  due  to  an  angular  spread  is 
skewed,  there  is  no  reason  to  expect  the  antisymmetric  properties  of  G(a9  =  0,v0)  to  be 


maintained  at  ae  >  0  .  As  seen  in  Fig.  1  and  (18),  the  characteristic  value  for  the  degradation  of 
gain  is  cr9*=1.  The  phase  velocity  for  peak  gain,  and  peak  absorption,  both  increase  with 
increasing  cr9  roughly  as  v0*  =  c9.  Note  that  the  general  features  of  Figs.  1  and  2  are  quite 
different  owing  to  the  different  forms  of  the  electron  distributions.  This  emphasizes  the  importance 
of  the  shape  of  the  electron  distribution  in  evaluating  gain  degradation  in  FELs,  and  the  need  for 
an  accurate,  flexible  theory  as  presented  here. 
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Fig.  2.  Intensity  and  contour  plot  of  ln(1+G(a0,vo))  with  j  =  1  and  a  =  0.  The  weak-field  gain 
degradation  in  this  low-current  FEL  is  due  to  an  electron  beam  angular  spread  with  a  normal 
distribution;  the  resultant  phase-velocity  distribution  is  the  exponential  distribution  function. 


VI.  High-Current  FELs 

The  integral  representation  of  the  optical  field  in  (15)  is  also  valid  for  high-current  FELs 
where  /  »  1.  in  this  case,  a(x)  acquires  a  non-linear  dependence  on  j  (  recall  expression  (7)  ) 
and  cannot  be  removed  from  the  integrand  of  (15).  To  proceed  analytically,  it  is  convenient  to 
remove  one  integral  from  (15)  by  taking  the  T-derivative  of  both  sides,  and  use  the  form  a  =a0eaz 


lor  the  optical  field.  Since  j  »  1,  assume  that  a  has  some  large  real  part,  even  though  the 
exponential  growth  may  be  somewhat  diminished  by  the  presence  of  cr  and  a0.  The  magnitude  a  0 
cancels  on  both  sides,  and  a  change  of  variables  gives  the  form 


The  upper  integration  limit  in  (19)  has  been  extended  to  infinity  because  the  integrand  containing 
the  factor  e~s  becomes  negligible  for  large  s . 

Eqn.  (19)  describes  several  properties  of  high-gain  FELs  without  integration.  If  cr,  v0,  and  a6 
all  -» 0  ,  then  a  has  the  same  roots  found  in  (7).  If  the  current  density  j  -» «>  so  that  a  real  part  of 
a  — *  oo,  then  we  obtain  the  same  limit,  since  a,  v0,  and  cre  all  appear  divided  by  a  in  (19).  Unlike 
the  low-current  FEL,  the  importance  of  beam  quality  in  a  high-current  FEL  depends  on  the  current 
density  j  .  This  feature  has  been  seen  in  FEL  experiments  and  simulations,  but  is  now  expressed 
analytically.  The  importance  of  beam  quality  can  be  made  more  quantitative  by  iterating  (19). 
Estimating  the  real  part  of  the  fast-growing  root  as  a*  =  ( j/2  )1/3V3/2,  the  integrand  of  (19)  is  only 
significantly  modified  when  a  *  =  ( j/2  )1/3  '13/2  or  when  cr9*  =  ( j/2  )v3^3/2.  In  the  high-current 
FEL,  the  characteristic  values  of  beam  quality,  a  *  and  cr9*  ,  are  not  equal,  and  increase  with  the 
current  density  j.  These  expressions  should  be  helpful  in  designing  high-gain  experiments  where 
there  is  a  trade-off  between  beam  quality  and  beam  current. 

Figure  3  shows  a  combined  intensity  and  contour  plot  of  In  (1+G  (a,v0))  for  moderately  high 
current  j  =  100  and  ct9  =  0.  The  points  at  the  peak  gain  In  (1+G )  =  4.3  are  indicated  by  white  on  the 
(a,v0)-surface;  black  indicates  zero  gain.  Contours  of  constant  gain,  In  (1+G)  =  2.0,  2.5,  3.0,  3.5, 
and  4.0  are  superimposed  on  the  intensity  plot.  For  the  high  quality  electron  beam,  small  cr,  gain 
is  confined  to  a  region  near  resonance,  but  extends  to  a  broader  range  in  v0  than  in  the  low- 
current  cases  of  Figs.  1  or  2.  This  agrees  with  the  discussion  below  (11),  and  gives  the  range  of 
optical  wavelengths  over  which  there  is  significant  gain  8v0  =  4 j  ~7  .  To  find  the  range  of 
wavelengths,  use  AX  fk ~  Av0  /2nN  about  the  resonant  wavelength  X  =  Xo(1+AT 2  )/2y2  .  The 
maximum  available  gain  decreases  significantly  as  a  ->  o  *  =  4.5  as  predicted  in  the  previous 
paragraph,  and  the  phase  velocity  for  peak  gain  roughly  follows  v0*  =  a  . 
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Fig.  3.  Intensity  and  contour  plot  of  In  (1+G  (a,v0))  with  =  100  and  cr9  =  0.  The  weak-field  gain 
degradation  in  this  FEL  is  due  to  an  electron  beam  with  a  normal  distribution  in  energy. 

Figure  4  shows  the  plot  of  ln(1+G(ae,v0))  with  c  =  0  so  that  the  gain  degradation  is  caused 
by  an  angular  spread  in  the  electron  beam.  The  contours  of  constant  gain  differ  from  Fig.  3 
because  of  the  new  shape  of  the  electron  distribution  function.  As  c9  increases,  there  is  a  slower 
decrease  in  gain  because  a9*  >  a*  as  found  above.  When  expressed  in  dimensionless  form,  an 
angular  spread  is  better  tolerated  in  an  FEL  than  is  an  energy  spread.  The  points  of  peak  gain 
increase  with  increasing  a9  similar  to  Fig.  3. 


Fig.  4.  Intensity  and  contour  plot  of  In  (1 +G  (a9,v0))  with  y  —  1 00  and  a  -  0.  The  weak-field  gain 
degradation  in  this  FEL  is  due  to  an  electron  beam  with  an  angular  spread  producing  an 
exponential  distribution  in  phase  velocities. 

Figure  5  shows  the  combined  intensity  and  contour  plot  of  ln(1+G(a,v0))  for  high  current 
j  =  io4  with  a0  =  O.  There  are  no  negative  gain  regions,  and  the  available  peak  gain  is  much 
larger  than  for  the  lower  current.  For  a  =  0.  the  position  of  peak  gain  is  essentially  at  resonance 
vo  =  0,  but  again  increases  roughly  as  v0*  =a  while  beam  quality  diminishes.  The  width  of  the 
gain  spectrum  at  a  =  0  is  wider  than  the  lower  current  case,  and  agrees  well  with  5v0  =  4 j  =12. 

The  contours  of  constant  gain,  ln(1+G)  =  14 . 24,  show  that  the  range  of  wavelengths  for  gain 

becomes  narrower  as  a  increases,  and  the  maximum  available  gain  decreases  significantly  as 

a  — >  a  *  =20. 
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Fig.  5.  Intensity  and  contour  plot  of  ln(1+C  (o,v0))  with  j  -  104  and  cr8-  0.  The  weak-field  gain 
degradation  in  this  high-gain  FEL  is  due  to  an  electron  beam  with  a  normal  distribution  in 
energy. 

Figure  6  shows  the  plot  of  ln(1+C(c9,v0))  for  high  current  j  =  104  with  cr  =  0.  Again,  the 

contours  of  constant  gain,  ln(1+G)  =  14 . 24,  are  distinct  from  Fig.  5  showing  the  importance  of 

the  electron  beam  distribution  function  even  at  high  gain.  As  c9  -» a9* ,  the  gain  decreases 
significantly,  but  again  the  angular  spread  is  seen  to  be  less  harmful  than  an  energy  spread. 
Unlike  Fig.  5,  the  position  of  peak  gain  stays  closer  to  resonance  as  o9  increases. 
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Fig.  6.  Intensity  and  contour  plot  of  In  (1+G  (o9,v0))  with  j  =  104  and  cr=  0.  The  weak-field  gain 
degradation  in  this  high-gain  FEL  is  due  to  an  electron  beam  with  an  angular  spread 
producing  an  exponential  distribution  in  phase  velocities. 
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ABSTRACT 

In  high-power  tree-electron  lasers,  the  electrons  can  become  trapped  in  deep 
potential  wells  formed  by  the  combined  optical  and  undulator  field  forces.  The 
trapped  current  oscillates  at  the  synchrotron  frequency,  and  can  drive  the  optical 
wave  at  sideband  frequencies  around  the  fundamental.  This  "trapped-particle 
instability"  can  occur  in  both  the  oscillator  and  amplifier  configurations. 

1.  Introduction 

Many  free-electron  laser  (FEL)  oscillators  [1],  and  FEL  amplifiers  [2]  are  designed  to 
produce  large  optical  power.  In  both  configurations,  a  relativistic,  high-current  electron  beam 
amplifies  a  copropagating  electromagnetic  wave  as  they  pass  through  a  periodic,  transverse 
magnetic  field  undulator.  The  trapped-particle  instability  starts  when  electrons  become  trapped  in 
deep  potential  wells  formed  by  the  combined  optical  and  undulator  fields.  Electrons  near  the 
bottom  of  the  well  oscillate  in  harmonic  orbits  at  the  synchrotron  frequency  causing  the  carrier 
wave  to  develop  sidebands.  The  "trapped-particle  instability"  was  first  predicted  using  a  single¬ 
mode  theory  describing  high-power,  low-gain  FEL  oscillators  with  a  tapered  undulator  [3];  the 
tapered  undulator  design  reaches  high  power  by  trapping  electrons  in  potential  wells  that  remain 
resonant  as  the  electrons  lose  energy  [4],  At  the  same  time,  multimode  simulations  of  short-pulse 
FEL  oscillators  observed  the  effects  of  the  instability  as  pulse  modulation  [5-7],  Recently,  the 
instability  has  been  observed  in  the  high-power  FEL  oscillator  at  LANL  [1],  and  possibly  in  the 
TRW/Stanford  oscillator  (8].  It  has  been  termed  the  sideband,  Kroll-Rosenbluth,  synchrotron, 
Raman,  and  the  trapped-particle  instability  [9-23].  Here,  the  multimode  sideband  theory  is 
reviewed  with  examples  from  short-pulse  FELs,  oscillators,  and  amplifiers.  Common  features  and 
differences  are  discussed. 
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2.  Multimode  Simulation  Theory 

The  electrons  injected  into  the  FEL  undulator  evolve  due  to  the  combined  optical  and 
undulator  electromagnetic  fields.  When  the  relativistic  electrons  ot  energy  Vo me2  (Yo»1)  are 
properly  injected  into  the  undulator,  the  transverse  electron  motion  is  periodic  with  amplitude 
)/27ty0  and  Phase  exp (ikQZ)  where  K  =  eB\0/2nmc2,  e  is  the  electron  charge  magnitude,  B  is 
the  rms  undulator  field  strength,  X0  =  2n/kQ  is  the  undulator  wavelength,  m  is  the  electron  mass, 
and  c  is  the  speed  of  light  in  vacuum.  The  average  speed  of  the  beam  along  the  Z-axis  is  p0c 
where  J30  =  1  ~('\+K2)i  We  follow  the  beam  evolution  with  the  dimensionless  time 
x  =  p 0ct/L  =  0-+1  along  the  undulator  length  L  =NXcw\hN  periods. 

The  light  wave  evolves  in  the  presence  of  the  electron  beam  according  to  the  transverse 
wave  equation.  The  optical  carrier  wave  with  frequency  c o  =  kc  =2ncfk  has  the  single-mode 
phase  exp[/(jfcZ-a>0]  and  a  complex  slowly-varying  coefficient  a{z)  =  aK(z)  +  ia,(z)  =  \a(z)  |  e ,<><z) 
evaluated  at  many  discrete  sites  z.  The  dimensionless  coordinate  z  is  the  ratio  Z/Nl  where  the 
"slippage  distance"  NX  is  defined  by  the  number  of  optical  wavelengths  that  pass  over  a  resonant 
electron  as  that  electron  traverses  the  undulator  length  L .  With  the  slowly-varying  amplitude  and 
phase  approximation,  and  the  coordinate  change  Z->Z+ct ,  the  wave  operator  reduces  to  a  single 
derivative  in  time.  The  dimensionless  optical  field  envelope  is  a(z)  =  4nNeKLE(z)lylmc2  where 
E(z)  is  the  complex  optical  electric  field,  and  y0 me2  is  the  resonant  electron  energy, 
Yo  =  k^+K2)/2k0.  Simulations  take  place  within  a  window  of  width  W  along  z  that  is  an  integral 
number  of  slippage  distances  long.  At  any  x,  electrons  at  coordinate  z  +  x  in  the  electron  beam 
overlap  the  light  at  coordinate  z  in  the  optical  wave.  The  electron  phase  relative  to  the  optical 
wave  and  undulator  fields  is  £  =  Ofc+/fc0)Z- cor ,  and  the  phase  velocity  is  v  =  £  =  L[(*+,fc0)p,  -k  ]. 


The  self-consistent  electron  and  optical  equat,.ns  [24]  are 

V(z+x)  =  [  1  -  V  (z+x)/2tuV  j  2  ( aR{z }cos(  C  (z+x)  )-a,(z)sin(  £  (z+x) ) )  (1) 

£(z+x)  =  v(z+x)  (2) 

aR(z)  =  -;(z+x)  <  (  1  -  v  {z+\)l2nN  ] 1,2  COS(  £  (z+x) )  >  (3) 

d,{z)  =y  (z+x)  <  [  1  -  v  (z+x)/2rtN  ] 1,2  sin(  £  (z+x) )  >  (4) 


where  the  dimensionless  current  density  is  j{z)~  8N(enKL)2p(z)/y$mc2,  p(z)  is  the  actual  particle 
density  at  site  z,  time  derivatives  are  with  respect  to  the  dimensionless  time  x,  and  <...>  is  an 
average  over  sample  electrons  at  site  z+x  in  the  electron  beam. 
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In  order  to  maintain  strong  coupling  between  the  electrons  and  light,  the  phase  velocity  v 
cannot  be  too  far  from  resonance,  v  =  0.  The  initial  resonance  condition  v  (x=0)  =  v0  is  taken  to  be 
at  the  maximum  of  the  weak-field  gain  spectrum.  In  both  the  amplifier  and  oscillator  cases,  we 
start  the  optical  field  at  low  values  a(x= 0)  =  aQ,  well  below  saturation.  The  dimensionless  current 
density  j  represents  low  single-pass  gain  when  j  =1 ,  and  high  gain  when  j  »1 .  Typically,  the  FEL 
oscillator  configuration  uses  j= 1,  while  the  FEL  amplifier  configuration  uses  y'»l.  Weak  optical 
fields  give  values  of  |a  |  <rr,  while  strong  fields  give  values  of  \a  \  >n.  For  small  current  and  weak 
fields,  maximum  gain  is  at  v0  =  2.6;  for  large  current  and  weak  fields,  maximum  gain  is  at  v0  =  0. 
The  electron  beam  at  each  z  has  a  uniform  distribution  of  initial  phases  £  (x=0)  over  a  2tc  range  so 
that  the  initial  optical  phase  is  arbitrary;  we  take  $(x=0)  =  0  at  each  z.  The  factors  [1-v/2nW]  in 
(1  )-(4)  are  close  to  unity  unless  the  electrons  lose  a  significant  fraction  of  their  energy  and 
become  less  relativistic.  Electron  motion  in  terms  of  v  corresponds  to  changes  in  the  electron 
energy  through  the  relation  v  =  4kN (  7-  Yo ) 1  To-  We  take  N  =  50  in  our  examples  as  a  common 
value  typical  of  both  amplifiers  and  oscillators.  The  final  optical  power  spectrum  P(k)  is  made 
more  relevant  by  expressing  k,  the  optical  wavenumber,  in  terms  of  the  corresponding  resonance 
parameter  v  {k );  similarly,  the  electron  distribution  function  /  (7)  is  expressed  in  terms  of  v  (y).  The 
power  spectrum  P(v(k))  and  the  electron  distribution  function  /(v(y))  are  more  physically 
meaningful  in  terms  of  their  affect  on  the  resonance  condition. 

The  single-mode  version  of  (l)-(4)  is  obtained  by  removing  all  z  dependence;  all  sites  are 
given  the  same  initial  values  of  C#,  v0,  and  a(z)  =  a0.  An  estimate  of  the  trapped-particle  motion  is 
obtained  by  considering  the  electrons  in  harmonic  orbits  near  the  stable  fixed-point  =  n/2.  With 
j  not  too  large,  and  and  small  energy  extraction  (N»  1),  \a  |  remains  approximately  constant 
during  the  synchrotron  oscillations  at  saturation.  Then,  the  motion  of  a  trapped  electron  is 
£(x)  =  C  +  ( vo  lvs)  sin(  v$x )  with  the  initial  position  (  £*,  v0).  The  synchrotron  or  trapped-particle 
oscillation  frequency  is  vs  =a^'2 .  When  the  trapped  electrons  oscillate  through  a  synchrotron 
cycle,  part  of  the  current  driving  the  optical  field  in  (3)  and  (4)  also  oscillates  through  one  cycle.  It 
is  the  oscillation  of  the  driving  phase  in  the  average  <...>  that  causes  the  trapped-particle 
frequency  to  be  imposed  on  the  optical  wave  as  it  slips  over  electrons.  The  sidebands  appear  at 
v0±  v5,  so  that  the  new  FEL  power  is  shifted  from  the  fundamental  wavelength  by  Ak/X  =  vs/2tdV. 
The  shift  has  a  simple  interpretation;  AX/A.  =  "the  number  of  synchrotron  oscillations"/  "the  number 
of  undulator  periods". 
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3.  The  Trapped-Particle  Instability  in  Short-Pulse  Oscillators 

An  FEL  oscillator  that  is  powered  by  an  RF  accelerator  injects  a  series  of  short  picosecond 
electron  pulses  into  the  undulator  while  the  optical  pulse  bounces  between  mirrors  separated  by  a 
distance  S  >L .  High-power  saturation  is  reached  after  many  passes,  and  the  FEL  continues  to 
work  for  an  additional  103  to  104  passes.  The  current  density  of  each  short  pulse  j(z)  is  taken  to 
be  parabolic  with  the  form  j{z)=j  (1-2 z 2/a2 )  for  | z  |  <  a2/^2  and  j (z )=0  for  | z  |  >  azH 2;  the 
length  az  is  normalized  by  the  slippage  distance  NX.  Typically,  RF  accelerators  produce  current 
densities  which  give  values  of  j  in  the  moderate  range  l-»  100  and  c,  =  1-»  30.  The  loss  on  each 
pass  (due  to  mirror  absorption  and  transmission)  is  described  by  e~nlQ  where  n  is  the  pass 
number.  Usually,  Q  is  from  2->  200. 

The  repetition  frequency  of  successive  electron  pulses  must  be  matched  to  the  bounce 
frequency  of  the  light  pulse,  2 Sc.  When  synchronized,  each  new  electron  pulse  arrives  at  the 
beginning  of  the  undulator  simultaneous  to  the  rebounding  optical  pulse.  The  "desynchronism," 
d  =  2AS/NX,  is  the  displacement  between  the  pulses  after  each  pass  when  the  mirrors  are 
separated  by  5  -  AS.  If  \d\  is  too  large,  the  electron  and  optical  pulses  do  not  overlap  for  a 
sufficient  number  of  passes  and  the  FEL  operates  below  threshold  coupling.  If  d  -  0,  exact 
synchronism,  the  FEL  is  also  below  threshold  (5-7,11-14,21]  due  to  an  effect  termed  "laser 
lethargy"  [12].  Because  of  slippage,  gain  is  preferentially  deposited  on  the  trailing  edge  of  the 
optical  pulse  causing  the  optical  pulse  centroid  to  travel  slower  than  c  in  vacuum;  therefore,  the 
optical  bounce  frequency  2Sc  is  overestimated.  To  compensate  for  the  "lethargic"  light,  the  path  5 
must  be  reduced  by  operating  at  d>0. 
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Fig.  1.  Driven  by  the  short  electron  pulse  j(z+ x),  the  trapped-particle  instability  has  modulated 
the  optical  pulse  envelope  a(z)  at  the  synchrotron  period  vs  =  2n.  The  total  optical  power 
P(n)  and  gain  G(n)  reach  steady-state  after  n=400  passes  through  the  oscillator.  G(v)  is 
the  weak-field  gain  spectrum,  included  for  reference.  The  final  optical  power  spectrum  P  (v) 
has  an  additional  sideband  at  v0  +  vs  due  to  the  instability.  The  final  electron  distribution 
/  (v)  is  broad  because  of  the  strong  fields. 

A  short-pulse  simulation  is  shown  in  Fig.  t .  Steady-state  is  reached  after  n  =  400  passes 
with  moderately  high  current  j  =  20,  an  intermediate  pulse  length  cr,  =  12,  N  =  50,  desynchronism 
d  =  0.04,  and  <2=3.  Large  output  coupling  is  used  (  low  <2  )  because  the  gain  is  far  above 
threshold.  The  top  right  figure  follows  the  power  P  (n)  averaged  over  the  whole  window.  Second 
on  the  right  is  the  net  gain  G  (n)  experienced  by  the  pulse  each  pass.  The  lower  left  picture  shows 
the  electron  pulse  shape  before  (x  =  0)  and  after  (x=  1)  the  slippage  of  one  unit  per  pass  in  a 
window  of  width  W  =  20.  The  upper  left  figures  show  the  final  optical  field  amplitude  \a(:)\  and 
the  optical  phase  profile  <j>(z)  plotted  along  z  where  -W/2<z  <W/2.  The  jumps  in  phase  of  2rc  are 
meaningless,  but  the  slope  of  the  phase  profile  indicates  a  local  change  in  the  laser  light 
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wavenumber  and  a  resonance  change  Av  (k )  =  -  $'(z )  away  from  v0.  The  grey  scale  on  $(z) 
shows  the  local  frequency  component  matching  the  gain  spectrum  G  (v)  on  the  right.  The  grey 
scale  superimposed  on  the  field  amplitude  shows  the  local  gain  experienced  by  each  part  of  a(z) 
on  the  last  pass  in  the  undulator.  This  grey  scale  matches  the  scale  in  G  ( n )  pictured  on  the  right. 
Third  on  the  right  is  the  final  electron  distribution  /  (v),  and  fifth  is  the  final  optical  power  spectrum 
P{v)  found  from  the  Fourier  transform  of  a(z).  The  fourth  picture  G(v)  is  the  weak-field  gain 
spectrum  for  the  current  density ;  =  20. 

The  simulation  has  reached  saturation  with  peak  fields  near  |  <2 1  =60.  Each  pass,  those 
electrons  near  the  center  of  the  pulse,  that  become  trapped,  encounter  fields  strong  enough  to 
cause  about  one  synchrotron  oscillation  since  ~  2n.  Normal  saturation  takes  place  early  and 
the  power  is  steady  for  about  50  passes;  then  the  power  increases  again  as  the  sideband  grows. 
Strong  optical  fields  and  the  resulting  trapped-particle  oscillations  must  occur  before  there  is 
significant  gain  at  the  sideband  frequency  vs.  The  extra  power  is  in  the  sidebands  and  the 
fundamental  remains  saturated  [17,23].  A  simple  calculation  shows  that  the  field  component  at 
the  sideband  frequency  grows  as  a(x)  =  as  exp  ((j/4vs) 1/2 1)  where  vs  =  \a  | 1/2  at  saturation,  and 
as  is  the  initial  sideband  field  at  vs.  In  this  simulation,  (j74vs)1/2  =  0.8  and  the  loss  rate  is 
-  (2  Q  )-1  =  -  0.2,  so  the  sideband  gain  is  above  threshold.  The  sideband  structure  is  clear  in 
I  <3  (z )  | ,  and  the  power  spectrum  P(v).  The  modulation  length  is  close  to  the  slippage  distance 
indicating  \s  =  2k,  and  the  sideband  in  P(v)  occurs  at  Av  =  2k  above  the  fundamental.  The  full- 
width  of  the  electron  distribution  is  given  by  the  height  of  the  closed-orbit  region  in  phase-space, 
4 1  a\v2. 

The  characteristics  of  the  trapped-particle  instability  in  short-pulse  FELs  are  briefly  outlined 
below  [21].  Most  characteristics  have  now  been  confirmed  by  experiments  [1 ,8,1 1]. 

1 .  At  small  d  >0,  the  FEL  usually  reaches  power  levels  large  enough  to  cause  the  trapped- 
particle  instability.  This  gives  a  broad,  possibly  chaotic,  optical  power  spectrum,  and  a 
broad  electron  distribution. 

2.  At  large  d,  the  steady-state  power  is  small  due  to  the  reduced  coupling,  and  the  trapped- 
particle  instability  is  less  likely  to  occur.  The  final  optical  power  spectrum  is  narrow  and  in  a 
single-mode;  the  final  electron  distribution  is  narrow  due  to  the  weak  optical  fields. 
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3.  When  d  is  in  the  intermediate  range,  we  have  often  observed  limit-cycle  behavior  in  the 
simulations  [5,6,13,21].  In  this  case,  the  pulse  continually  changes  shape  while  the 
trapped-particle  instability  creates  new  subpulses. 

4.  Increasing  the  current  density  j  or  the  resonator  Q  increases  the  steady-state  power,  the 
synchrotron  frequency,  and  the  sideband  gain.  The  addition  of  sideband  power  is 
cumulative,  since  the  presence  of  a  strong  sideband  again  increases  the  steady-state 
power. 

5.  If  taper  is  introduce  into  the  undulator  design,  the  synchrotron  frequency  is  only  slightly 
modified,  and  the  sideband  gain  is  reduced  [13,22,23]. 

6.  When  the  pulse  length  at  is  near  unity,  the  short-pulse  effects  can  be  dramatic  [11].  Since 
the  sideband  modulation  first  appears  near  the  slippage  distance,  a  short  pulse  may  not  be 
long  enough  to  support  modulation  at  the  synchrotron  frequency.  The  optical  pulse  can  be 
made  significantly  shorter  or  longer  than  the  picosecond  electron  pulse  by  adjusting  d. 

4.  Simulations  in  a  "Wrapped-Window" 

If  the  electron  pulse  is  not  short  (<rt>>1)  it  is  prudent  to  simulate  the  FEL  by  sampling  a 
smaller  window  W <<cr2  with  periodic  boundary  conditions  such  that  £  ( z-Wl2 )  =  £  (z+W/2)\  these 
can  be  called  "wrapped-window"  simulations  [17,22,23].  Using  a  number  of  sites  Nw  in  a  window 
of  width  w  we  follow  a  restricted  number  of  modes  v,  =  v0-  {2n/  W){  l-Nw/2 )  where 

/  =  0, 1 , 2 . Nw-1 ;  the  mode  spacing  is  Av  =  2rtW .  The  desynchronism  d  and  the  pulse  length 

cr2  don’t  enter  the  problem,  and  the  current  density  j  (z  )=j  is  constant  along  the  window. 

The  FEL  simulation  in  Fig.  2  uses  j  =20,  2=3,  and  N  =50  (as  in  Fig.  1),  but  for  a  long 
pulse  a2>>W=4.  The  individual  pictures  are  the  same  as  in  Fig.  1  except  for  the  reference  to 
pulse  shape.  Without  some  noise  source,  no  power  would  develop  at  frequencies  other  than  the 
fundamental  and  every  site  z  would  evolve  identically.  Specific  sources  of  noise  can  vary  from 
one  experiment  to  the  next.  In  Fig.  2,  the  initial  electron  phases  are  uniformly  spread  over  a  2 n 
range,  but  with  an  additional  random  phase  of  zero  mean  and  standard  deviation  8£  =  l.OxiCf6. 
Electron  shot  noise  is  a  typical  source  of  this  kind  of  noise.  In  the  pulse  simulations  no  random 
noise  is  introduced,  since  the  spectral  features  of  j(z)  are  sufficient.  After  n  =400  passes,  the 
steady-state  optical  fields  peak  at  |  a  |  =  60  just  as  in  Fig.  1 ,  and  the  final  power  spectrum  P  (v)  is 
similar  to  Fig.  1 ,  but  without  the  short  pulse  features. 


Fig.  2.  The  "wrapped-window"  simulation  gives  a  clear  example  of  the  trapped-particle  instability 
in  the  FEL  oscillator  without  the  complicating  features  of  short  pulses. 

Wrapped-window  simulations  have  many  of  the  same  characteristics  as  for  pulses,  but 
without  the  effects  of  desynchronism.  As  a,-»  °°  the  wrapped-window  approach  is  the  only  viable 
solution  to  the  FEL  multimode  simulation. 

1 .  Increases  in  j  or  Q  make  the  trapped-particle  instability  more  severe.  It  is  unusual  to  find 
the  periodic,  even  modulation  shown  in  Fig.  2. 

2.  When  there  is  taper  in  the  undulator  design,  the  instability  is  less  prominent  for  the  same 
reasons  as  stated  in  the  pulse  case. 
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3.  Over  many  hundreds  of  passes,  the  stored  optical  wave  "sees"  many  synchrotron 
oscillations,  so  that  any  sideband  gain  above  threshold  gives  large  growth  from  a  small 
amount  of  noise.  The  resulting  steady-state  features  are  therefore  not  affected  by  the 
details  of  the  noise  source  employed. 

5.  The  Trapped-Particle  Instability  in  High-Current  FEL  Amplifiers 

In  the  high-current  FEL  amplifier,  y>>1,  large  optical  power  can  be  produced  in  a  single 
pass  through  the  undulator.  An  induction  linac  accelerator  and  a  long  undulator  can  result  in 
values  j  =  104— »105  [2J.  The  electron  pulses  are  long,  so  that  the  wrapped-window 

simulation  method  is  essential.  The  FEL  growth  rates  are  so  large  that  electrons  become  trapped 
early  in  the  undulator  and  begin  executing  synchrotron  oscillations.  There  can  be  from  several  to 
many  tens  of  synchrotron  oscillations  along  the  undulator  so  that  the  optical  field  experiences  far 
fewer  synchrotron  cycles  than  in  the  oscillator  case.  However,  even  the  limited  number  of 
synchrotron  cycles  can  result  in  significant  sideband  gain  owing  to  the  large  current  density  j . 

The  wrapped-window  simulation  in  Fig.  3  follows  the  power  and  gain  from  t  =  0->1  for  an 
FEL  with  high  current  j  =8xl04,  an  initial  field  a0  =  20  starting  at  vo  =  0,  and  N  =  50.  No  electron 
phase  noise  is  present,  but  a  small  initial  sideband  field  as  =  0.01  is  introduced  at  vs  =  207t. 
Spontaneous  emission  or  electron  shot  noise  from  the  trapped  electrons  can  lead  to  a  contribution 
of  this  size  because  j  is  large.  The  power  />(t)  and  gain  G(t)  are  seen  to  oscillate  with  the 
synchrotron  frequency  vs  =  20rc  and  impose  a  strong  modulation  on  the  light  wave  envelope  a  (z ). 
The  final  power  spectrum  P(v)  shows  that  the  sideband  has  grown  to  almost  equal  the 
fundamental  power. 


- 10- 


*909610*  FEL  Optical  Modes 

j=80000  aQ= 20  N=50 

as^.Oi  ks=S2.B 


Fig.  3.  The  trapped-particle  instability  can  create  significant  sideband  power  during  one  pass  in 
the  FEL  amplifier.  The  high-currents  used  in  amplifiers  give  much  stronger  optical  fields  at 
early  stages  of  the  undulator  that  can  cause  several  synchrotron  oscillations,  and  impose  a 
corresponding  modulation  onto  the  optical  wave  envelope.  In  this  example,  the  sideband  is 
spaced  at  vs  =  20 jc. 

The  trapped-particle  instability  in  FEL  amplifiers  differs  in  many  ways  from  the  oscillator 
case.  Sources  of  noise  are  much  more  important,  and  there  is  no  resonator  Q  to  consider. 

1.  Large  current  j  makes  the  trapped-particle  instability  more  likely.  Using  the  sideband 
growth  rate,  we  can  evaluate  the  current  j  needed  to  make  the  sideband  field  as  equal  to 
the  fundamental,  a0  =  v|.  This  gives 

/  =  I6vjln2(vj/Va7)  . 
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The  current  j'  has  a  slow  logarithmic  dependence  on  the  initial  sideband  field  as  so  that  an 
accurate  estimate  is  not  too  crucial;  we  take  as  =  0.01 .  We  also  estimate  that  there  must  be 
at  least  a  few  synchrotron  oscillations,  say  vs  =  l6rt,  before  the  optical  wave  could 
accurately  determine  the  synchrotron  frequency.  This  gives  a  characteristic  current  density 
j'  =3xl04  where  sidebands  could  be  expected  to  significantly  alter  the  final  FEL  amplifier 
spectrum. 

2.  In  the  tapered  undulator  case,  less  current  is  trapped,  and  the  synchrotron  oscillations  have 
a  smaller  amplitude  so  that  the  trapped-particle  instability  is  observed  to  be  less  severe,  as 
in  FEL  oscillators. 

3.  Input  noise  at  the  sideband  frequency  is  important  to  the  development  of  significant 
sideband  power  in  FEL  amplifiers.  Shot  noise  and  spontaneous  emission  cannot  be 
eliminated,  and  may  be  a  significant  contribution  at  large  current  densities. 

4.  For  the  large  currents  j>j‘ ,  the  FEL  does  not  reach  steady-state  operation  even  in  strong 
optical  fields  |  a  \ .  The  power  continues  to  increase  and  so  does  the  synchrotron  frequency. 
This  continual  change  in  the  synchrotron  frequency  may  play  a  useful  role  in  suppressing 
the  growth  at  any  particular  sideband  frequency. 

In  conclusion,  we  find  that  the  trapped-particle  instability  is  more  difficult  to  suppress  in  the 
FEL  oscillator  than  in  the  FEL  amplifier.  The  large  number  of  synchrotron  oscillations  experienced 
by  the  light  stored  in  the  oscillator  make  the  effective  gain  over  for  many  passes  large.  But,  while 
the  trapped-particle  instability  has  already  been  observed  in  the  FEL  oscillator,  it  has  not  yet  been 
found  in  the  FEL  amplifier.  The  many  differences  between  FEL  oscillator  and  amplifier 
simulations  implies  that  we  should  not  use  the  oscillator  case  as  a  direct  proof  that  there  will  be  an 
equally  prominent  instability  in  real  FEL  amplifiers. 
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